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In this paper we propose a stochastic broadcast ;r-calculus which can be used to model server-client 
based systems where synchronization is always governed by only one participant. Therefore, there is 
no need to determine the joint synchronization rates. We also take immediate transitions into account 
which is useful to model behaviors with no impact on the temporal properties of a system. Since 
immediate transitions may introduce non-determinism, we will show how these non-determinism 
can be resolved, and as result a valid CTMC will be obtained finally. Also some practical examples 
are given to show the application of this calculus. 

1 Introduction 

Process algebras such as CCS |[T91 . CSP IfTTl . and AC? 121 have been successfully used to model and 
analyze concurrent systems. The system behavior of these classical process algebras is usually given 
by labeled transition systems (LTS) which have proved to be a convenient framework for analyzing 
qualitative properties of large complex system. As these models are only concerned about functional 
aspects of concurrent systems, process algebras have been extended with stochastic variables in order 
to model performance-oriented systems in recent years. Such examples include TIFF [11], FEFA |[T5l . 
EMPA m, stochastic ;r-calculus llQl, IMC 1131, StoKlaim HI, and Stochastic Ambient Calculus ED. 
The semantics of these models are given by a variant of LTS, Continuous Time Markov Chain (CTMC), 
which can be used to analyze quantitative properties directly. Each transition in a CTMC is associated 
with an exponentially distributed random variable which specifies the duration of this transition. The 
underlying CTMC captures the necessary information for both functional verification and performance 
evaluation. 

Synchronization in stochastic scenarios have been addressed in |[T5l[T3l [TI using different techniques. 
In this paper we develop a stochastic broadcast ;r-calculus aiming at modeling server-client based sys- 
tems which are used widely in practice. In such systems synchronization are always governed by one 
participant, so there is no need to determine synchronization rates like others. In our calculus only out- 
puts are associated with rates and their durations are exponentially distributed while inputs are always 
passive. We all know that the nondeterministic choices among outputs can be resolved by race condi- 
tions probabilistically. Similarly, to resolve nondeterministic choices among inputs, we let each input be 
associated with a weight as usual and the probability of an input receiving a message is determined by its 
weight and the total weight of all current inputs. In addition the communication in our calculus is based 
on broadcast, that is, when one component outputs a message, it will be received by all the recipients 
instead of only one of them. Such scenarios can be found in practice very often. For example consid- 
ering the checking out in a supermai^ket, the arrivals of customers can be assumed to be exponentially 
distributed. When a customer comes to the counters, he/she will choose different counters according to 
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the lengths of their queues, the longer the queue the less likely it will be chosen, meanwhile when the 
customer is checking out, not only the counter knows it but also other departments will know and react 
accordingly such as financial, purchasing and so on. 

To enhance the expressiveness of our calculus we also take immediate actions into account. The 
immediate actions will happen instantaneously and have been studied in mdlKMl. They are useful to 
describe certain management and control activities which have no impact on the temporal behavior of 
a system. Since immediate action takes no time to execute, so race condition does not apply here. In- 
stead we assign each immediate action a weight to resolve nondeterministic choices between immediate 
transitions which is similar as inputs. For instance activities such as "when the buffer of a server is full, 
the coming clients will be transferred to another server instantaneously" can only be modeled by using 
immediate actions. In this paper we give several classical models from performance analysis which can 
be modeled in a compositional way by making use of immediate transitions. Accordingly, we will call 
the non-immediate transitions (resp. actions) Markovian transitions (resp. actions) in the sequel. 

Usually the problem of immediate action is that the existence of an underlying CTMC can no longer 
be guaranteed. In this paper, we solve this in two ways. As usual immediate transitions take no time 
and should have priority over Markovian transitions, so when an immediate transition is available, it 
will block the executions of Markovian transitions. We divide the whole process space into two sets: 
Immediate Processes (IP) and Markovian Processes (MP). IP only contains processes where at least one 
immediate transition is available and MP contains processes where no immediate transition is available. 
Since immediate transitions can exempt the execution of Markovian transitions, we can say that states in 
IP can only perform immediate transitions. All states in a CTMC will belong to MP. To calculate the 
rate from P to P' in a CTMC, we accumulate the rates of all the possible transitions from P to P' where 
transitions might be via states in IP. Sometimes it is possible for a process reaching a state which and 
all its derivations belong to IP. In this case, no time is allowed to elapse and the process is said to be 
absorbing. We use a special state Stuck to denote such situation and show how a CTMC can be obtained 
even with the existence of immediate actions. 

Similar with the existing calculi whose semantics are given by LTS, we also give the LTS for our 
calculus. Differently, each Markovian transition in our LTS is labeled by a rate instead of an action. 

For example, a typical transition looks like P — '—^ AP where A denotes that the execution time of this 
transition is exponentially distributed with rate A and AP is a distribution over pairs of action and process 

(a, Q). Intuitively, if P — ^^ AP, that means that P will leave its original state with rate A (sojourn time of 
P is exponentially distributed with rate A) and get to Q via action a with probability p if the probability 
of {a, Q) in AP is equal to p. By defining an LTS in this way, the correspondent CTMC can be obtained 
in a natural way. It is worth mentioning that our framework could also be used as an alternative general 
way to specify the LTS as rate-base transition systems ^. Without relying on different techniques, for 
example multi relations, proved transition systems and unique rate names used in PEPA, stochastic n- 
calculus, and StoKlaim respectively, we can have a uniform way to define the underlying models for 
these stochastic calculi. 

The paper is organized as follows: the syntax of our calculus is presented in the next section and 
in Section 3 we give the Labeled Transition System. In Section 4 we illustrate the use of immediate 
transitions by giving a few examples. We show how to get the underlying CTMC even with existence of 
immediate transitions in Section 5. Finally, we end by concluding and describing the future work. 
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2 Syntax 

Before introducing our calculus, we first give the following general definition of probability space. A 
probability space is a triplet P = (Q, F, 77) where O is a set, f is a collection of subsets of O that includes 
Q. and is closed under complement and countable union, and 77 : F ^ [0, 1] is a probability distribution 
function such that t]{Q.) = 1 and for any collection {C,}, of at most countably many pairwise disjoint 
elements of F, rjiUiCj) = YjiViCt)- A probability space {D.,F,r]) is discrete if fl is countable and F = 2", 
and hence abbreviated as (fl,?/). Given probability spaces [Pj = (fl,-,?7,)};e/ and weights w, > for each 
/ such that Yjiei'^i = 1' the convex combination Yjiei^fPi is defined as the probability space (fl,?/) such 
that Q - U/e/^r ^rid for each set Y QQ., r]{Y) = 2,g/W;77,(y nQ,)- Usually, we use [pi : Pi]iei to denote a 
probability space P = ({P/l/e/,?/) such that ?7({P;}) = p,, here / is a countable index set. Dirac probability 
space {I : P) will be written as P directly in the sequel. If YjieiPi ^ 1 then we call it a sub probability 
space. We also use P{Pi) - Pi to denote the probability of P, in P. The summation and pai^allel between 
two sub probability spaces can be defined in a natural way as follows: 

Pl+P2^ {pi +P2 : P I Pl(P) =PlAp2{P)^P2APmi)+P2(^2) < U, 
Pi II P2 - {Pl XP2 : Pi II P2 I PliPl)=Pl Ap2iP2)=P2]. 

Note in the above Pi{Q.i) + p2i0.2) < 1 is used to guarantee that P\ + P2 is still a valid sub probability 
space. 

We presuppose a countable set 'if of constants and a countable set Y of variables ranged over by 
a,b,c... and x,y,z--- respectively such that "^ n ^ = 0. n, m, / . . . e "^ U ^ are called names. The syntax 
of processes is given as follows where A e l.>o is the exponential rate and w e ]R>o is the weight of the 
input action. When the rate of an output is infinite, it is an immediate action which takes no time for 
it to be performed. We use cx3„. to denote an infinite rate with weight w. In the following A is used to 
denote either A or 00^. and Aq ranges over exponential rates as well as 0, that is, Aq e R>o. It is obvious 
that every output action must be prefixed by an exponential rate and every input action has a specified 
weight, if the rate of an output is infinite then it will be assigned with a weight instead. We assume that 
there is a countable set of constants, ranged over by A, which are used to denote processes. By giving an 

def 

equation such that A = P we say that constant A will behave as P, here A is required to be guarded in P, 
i.e. every constant appearing in P has to be prefixed by fKct. We only consider closed processes here and 
use P,Q-- to range over closed processes J^. 

Act ::= n(x,w) \ n{m,A) 

P, e ::- I kct.P \vaP\P+Q\[n^ m]P, 2 | P || 2 | A 

A substitution {a/x] can be applied to a process or process distribution. When applied to a process 
distribution, it means applying this substitution to each process with probability greater than in it. 
The set of free names and bound names in P, denoted by fh{P) and bn{P) respectively, ai^e defined as 
expected and n{P) -fn{P) Ubn{P) denotes the set of all the names in P. Structural congruence, =, is 
the least equivalence relation and congruence closed by the rules in Tableland a-conversion. = is also 
extended to network distributions as usual. 

3 Semantics 

The actions of processes ^, ranged by a,/3 • • • , are defined by 

a ::= a{x) \ vba{b) \ r 
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Table 1: Structural Congruence 



vavbP = vbvaP P + Q = Q + P P\\Q = Q\\P {a = a\P,Q = P 
[a = b]P,Q=Qa^b (vaP) \\ Q = va{P \\ Q) a ifn{Q) 



Table 2: Function y evaluating weight of input on a given channel 



y{a,0) = 

Iw Act - a(x, w) for some x 
otherwise 

(O a = b 

y{a,vbP) = { 

\y{a,P) otherwise 

y{a,P + Q) = y{a,P) + y{a,Q) 

iyia,P) b = c 
y{a,[b^c]PQ) - r ' 

[yia,Q) bi=c 

y{a,P\\Q) = y(a,P) + y(a,Q) 

def 

y{a,A) = y(a,P) A = P 



Here B is a set of constants, when b eb, bis bounded, otherwise it is free. The functions ^,Z7«, and n can 
be lifted from processes to actions as usual. 

To evaluate the total weight of inputs on a given channel in a process, we define function y : "^ x =^ — > 
R>o as Table |2] 

Similarly, we also give the function fi:'^x^—^ R>o to evaluate the total rate of outputs on a given 
channel in a process which is defined in the Table |3] To evaluate the weight of outputs with infinity rates, 
we define function //oo :^x ^ —^ R>o which is the same as /i except that: 

\w A = oo„ Aa - b 
IUco{a,b{ni,A).P) = < 

otherwise 



In addition, //(P) = YjaeVl^i^^P) ^^id jicaiP) - TjaeVf^ooia,P) are used to evaluate the total rate and 
total weight associated with infinite rates of outputs in a process. 

We define process distribution, ranged over by P,Q. . ., as a probability space where Q. = J^. Sim- 
ilarly, process action distribution can be defined as a probability space where Q. = Jix ^. We use 
AP,AQ... to range over process action distributions. The set of all the actions in AP is defined by 
J?[(AP) = {a I 3P.AV{a, P) > 0} while the corresponding sub process distribution of a in AP is denoted 

by APia) = {p : P \ A¥{a,P) = p > 0}. We will write AP as (a,P) if J?l(AP) = {a] where P = AP(a). In 
addition, we use AP(P) = ^ AF{a,P) to denote the total probabihty of P in AP. 

We lift new operator to process action distributions in ([T]). If the channel is restricted, then the 
broadcast action will change to r; if the message is restricted, then the broadcast action will be updated 
accordingly; otherwise the broadcast will stay unchanged while the new operator will be put on the result 
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Table 3: Function ju evaluating rate on a channel of a process 



Ka,0) = 

iA A = AAa = . 
otherwise 

lu{a,n(x,w).P) = 



/u{a,b{m,A).P) 



/j.{a, vbP) 



JO a = b 

\^{a,P) otherwise 
lu{a,P + Q) = n(a,P)+ii{a,Q) 



IJia,[b = c]PQ) = 



\fi{a,P) b-c 
\n{a,Q) b + c 
^i{a,P\\Q) = Ka,P)+l^{a,Q) 

H{a,A) = Ka,P) At^P 



process. 

vaAF ={p : (r, vaP) | AF{vba{b),P) = p] 

U {p : {vab{a),P) \ AF{b{a),P) =pAai=b} (1) 

U{p : {a,vaP) \ AF{a,P) = pAa ^fn{a)} 

The semantics of our calculus is shown in Table U] where / and / are finite index sets. We use 

-— '*- to denote *- or — ^ , where s- is the Markovian Transition with rate A, and - - ^^ is 

the Immediate Transition with rate infinity and weight w. A transition with rate *- is called a 

passive transition (TSl . All the transitions have the form P-— -^AP which means that P will evolve 
into process Q by performing action a with probability p if A¥{a, Q) = p. In addition, when it is a 
Markovian transition with rate A, it means that P will leave to other states with rate A or the duration 
of the transition is exponentially distributed with rate A. It is not hard to see from the semantics that 
for a Markovian transition, all the actions in the resulting distribution J?l(AP) are either outputs or t 
actions, while for the passive transitions, J?l(AP) only contains an input action, therefore can be written 
as (a{x),¥) where P = A¥(a{x)). Rule (REC) means that process a(x,w).P can receive a message on 
channel a and then evolve into P with probability 1. Similarly, in rule (mBRD) a{b,A).P will leave to 
other states with rate A and evolve into P by broadcasting the message b on channel a with probability 
1, this is a Markovian transition. If the rate of output is infinite, it should be performed instantly. This 
is called immediate transition which is shown by (IBRD). The weight associated with the infinite rate is 
used to resolve nondeterministic choices as in input actions. Rule (RES) only applies to Markovian and 
immediate transitions, since A>0 can guarantee that the transition is not passive. The new operator on 
process action distribution is defined by ([T|). By definition of y in Table [2l if y{a,P) - which means 
P is not ready to receive messages on channel a, in this case P will ignore all the messages broadcasted 
on channel a. This results in rule (LOS). Every input action is associated with a weight which can be 
used to resolve nondeterministic choices among diff"erent input actions probabilistically. For example 
after receiving a message b on channel a, a{x,2).P\ + a{x, \).P2 will evolve into P\{blx] with probability 
2^ and P2{blx] with probability ■:^. This is shown in (SUMl). (R\R1) is straightforward since our 
calculus is based on broadcast. Two parallelized processes will evolve together after receiving a message 
on a certain channel. Intuitively, when we put processes Pi and P2 together, the compositional process 
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P will leave to other states with rate Ai + A2 if the rates of Pi and P2 for leaving their original states 
are A\ and A2 respectively. Whether Pi or P? will be executed first depends on the race condition, that 
is, Pi will be executed before P2 with probability j-pj; and the probability for the other case is ^ '^^ . 

This is captured by rules (SUM2) and (PAR2) when -— -^ = ^ . In (PAR2) when -— -*. = ^ we 

also need to consider all the possible synchronization between P and Q. For example, if P can evolve 
into a sub process distribution AF{vbiai{bi)) after action vbjaiibi) and Q will evolve into Qi{bi/x] after 
receiving bi on channel a,, then P \\ Q will evolve into sub process distribution AF{vbiai{bi)) \\ Q,{Z7;/x} 
by performing action vbjai{bi) after leaving from the original state {vbi C\fn{Q) = 0). Since P and Q may 
have several outputs available at the same time, we need to list all the possible synchronization and then 
add all the resulting sub process action distributions to form the final result. The following example is to 
show how (PAR2) works. 

Example 1 Suppose we have two processes: P = n(y,3) \\ {m{x,2).P\ + m{x,A).P2) and 

3 
Q = m{z,2) II n{x,l).Qi. By the semantics, P ^{1 : {n{y),m{x,2).P\ + m(x,4).P2)} = {n{y),f) and 

2 
Q ^ {1 : {m{z),m{x, l).Q\)] = (m(z),Q). When we put the two processes in parallel, we have to con- 
sider all possible synchronization between them. P can broadcast y on channel n and Q can broadcast z 
on channel m, in the meanwhile P can receive a message on channel m and Q can receive a message on 
channel n, formally, 

U: {m{x)My,^)\\Pi)\_. . ,^. 
P ^ { 9 ) = (m(x),P ), 

(f: {m{x),n{y,^)\\P2)] 

e^U{l : (n{x)Mz,2) II Gi)) = («W,Q')- 

In P II Q, either P or Q will broadcast a message first, and the non-determinism is resolved probabilis- 
tically by race condition, i.e. P || 2 will perform n{y) first with probability t and the probability ofm{z) 
being executed first is j. When n{y) is performed, Q will receive it and evolve into '^{ylx}. Similarly, 
when m{z) is executed, P will evolve into P'{z/x} accordingly. So 

P II Q^^\^{n{y),¥ II Q'{3;/4) + ^ x(m(z),Q II P'{z/x)). 

Rules (SUM3) and (PARS) are similar with rules (SUM2) and (PAR2), but they only apply to pro- 
cesses where Q can only have a passive transition (with label 0), this is guaranteed by n{Q) = 0. (SUM3) 
and (PAR3) cannot be omitted since in (SUM2) and (PAR2) both P and Q are required to have non- 
passive transition, while in (SUM3) and (PAR3) only one of them has non-passive transition. The argu- 
ments for these rules when -— «» - — ^ are similar. Rules (CON) and (STR) are standard and need no 
more comments. 

From the syntax and semantics we know that the nondeterministic choices among Markovian tran- 
sitions can be resolved by a race condition while both the nondeterministic choices among immediate 
outputs and inputs can be resolved based on their weights. But still there might be nondeterminis- 
tic choices during the evolution of a process, such as nondeterminism between passive transitions and 
Markovian transitions and nondeterminism between immediate transitions and Markovian transitions. 
These nondeterminism can be resolved easily since we assume that immediate transitions can preempt 
other transitions while passive transitions should not be considered when talking about the underlying 
CTMC of a process. We will talk about this in details in Section [5] 

In this section we will not discuss immediate transitions, we leave it to the next section. The following 
simple example is to show how to get a CTMC from a process without immediate actions and we often 
omit the tail process 0. 
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Table 4: Inference Rules(^^ denotes either — > or -->) 



(REC) 



(mBRD) 



a{x,w).P—^{l :(fl(x),P)} 



a{b,A).P—^{l -.(aib^P)] 



(iBRD) 



(RES) 



a{b,oo,,}.P--^[l:{a{b),P)] vaP--^ 
,xifn{P) and y{a,P) = (LOS) 



P^^[l:{a{x),P)] 

P— ^(a(x),P) Q-^{a{x),q 
P + e ^^ (fl(x), rr^^i^O^P + ^^^^ 



,y{a,P + Q) + Q{S\5M\) 



Q 



P + Q 



''l+'*2 ,!| 



y(a,P + Q)'^ ^ y{a.P + Q)^' 

^^ (SUM2) 



/1 1+^2 



-^2 
^1+^2'' 



-^AMe)=o)v(- 



A;Uoo(e) = 0) 



(SUMS) 



p + e- 



P— ^(a(jt:),P) 2— ^(a(jt;),« 
^lie-^(flW,P||Q) 



(PARI) 



Q- 



U bj nfniP) = V vbjajibj) E 



U bi r\fii{Q) = V vbiOiibi) e 



».^— ^(«,W,P,) 



e — ^(fl,«,Q,) 



^ '*' ( + (vS,-fl;(^), AP(yS,-fl,0,-» II Q.{bi/x]) + (t, AP(t) || Q)) ^ 



(PAR2) 



^IIG 



,l,+,l2 



/il + /I2 ie/ 

+ — ^( + (vbjaj{bj),¥j{bj/x] II AQ(yV;<^-») + (r,P II AQ(t))) 
/il +/I2 ;£■/ 



''p-.^AP VvS,fl,<Z7,)ey[(AP).e— ^(a,(-*^),Q,) U^nX0 = 0'' 

16/ '£/ 



A;u(e)=o)v(- 



Ayu.o(e) = 0) 



P II Q -— ^ + (v/7,fl,<^/), AP(v^fl,-0/» II Q,{^,7-^}) + (r, AP(t) II 

i€l 



(PAR3) 



-io 



-io 



clef 

,A d p (CON) 



P^Q 



■h 



(STR) 
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r 



p 

^ y'Ve i/4\ ^ x^-"^/^ ^ v"""*""^ *^ 

Pn Pi2 P21 P22 Pii P12 Pi\ Pii 

(a) (b) 

Figure 1: A Simple CTMC 

Example 2 Given a process P = a{b\,2) || a{b2,6) || (a(x, l).Pi + fl:(x,2)./'2). ^''^'^ P will broadcast a 
message (b\ or b2) on channel a with an exponential delay 8 and then evolve into a process action 
distribution AP. By semantics in Tabled we know 

AP = {-i- : {a{bi},Pn),2- {a{bi),Pn),--. {a{b2),P2i),^: {a{b2),P22)} 

where 

Pn = a{b2,6) II Piibi/x] Pn = «<^2,6) II ^2(^1/4 
P21 - a{bu2) II Pi{b2/x] P22 = a{bu2) \\ ^2(^2/4 

This is displayed in Fig. [Tfaj and the correspondent CTMC is shown in Fig. \l\b), here we use dot lines 
to denote probabilistic choices and omit actions of the transitions. 

In the above example we briefly illustrated how to get a CTMC from a process. Now we are going to 
give the general construction by which we can get the coiTcspondent CTMC from a process. Use M{P, Q) 
to denote the rate from P to 2 in a CTMC, and define Der(P) as the smallest set of processes satisfying: 

i) P € Der(P); ii) P2 e Der(P) iff there exists Pi e Der(P) such that Pi ^^ AP with AP(P2) > 0. So 
Der(P) is the set of all the processes which are reachable from P with positive probability via arbitrary 
steps. For each two processes P\,P2 £ Der(P), the rate from Pi to P2 is equal to A x AP(P2), that is 

M{P\,P2) = ix AP(P2) such that Pi — ^AP, otherwise M(Pi,P2) = 0. 

4 Immediate Transitions 

In this section, we will give a few examples and show how can we benefit from immediate transitions. 

First we consider a model called Closed Queueing Networks (CQN) ||7l from performance analysis. 
A queueing network is a collection of servers. Customers must proceed from one server to another 
in order to satisfy their service requirements. The queueing network is closed if neither arrivals nor 
departures of customers are permitted; instead the number of customers in the network are fixed at all 
times. We use M to denote the number of servers in the network and N for the number of customers 
circulating around. The service time for a customer at server / is exponentially distributed with rate Ai 
and the probability a customer will proceed to the j-th server after completing a service request at server 
/ is equal to pij for /, 7 = 1,2,..., M. 

Example 3 Suppose we are given a CQN with 5 servers and 15 customers shown in Fig. ^where the 
numbers in the rectangles denote the length of the queue of each server as well as their indexes, the 
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Qi 



Go 



23 



Figure 2: A CQN with 5 Servers 



Figure 3: Open Queue Network with Blocking 



numbers on the edges denote the transition probabilities and the numbers in the circles denote the rate 
of each service time. 

The communication in the CQN is "point-to-point" in the sense that each leaving customer will 
arrive at only one server, while in our calculus the communication is based on broadcast. But with 
immediate transitions we can model such "point-to-point" communication as follows. Here we assume 
that the weight of each input is 1 by default, that is, Ci{x) is equal to Ci{x, 1). SQi{wi) denotes server i 
with Wi customers waiting for service in its queue. In Table\5\ CQN{ 1 , 2, 3, 4, 5) denotes the system where 
the i-th parameter is the length of the i-th queue. 







Table 5: 


Model of Closed Queueing 


Network 








Rev(/ 


j,Wj) = Ci(x).[x = w]{c{r,oop_^).(SQj(Wj + 1) | 


SQ,(w- 


-1)) + c(x).SQ;(h 


;■)) 


we 


[1,15] 


SQj(wj) = X 


Re\/{i,j,wj) + 


(K=0]0,c,<w 


^^j)) 














l<'<5,P/j 


>o 
















CQN(1,2,3,4,5) = 


vcdiii SQ,(0) 

















Each SQj(Wj) contains two parts: receivers denoted by ^i<i<5nj>of^SV{i,j,Wj) listen on the chan- 
nels of their predecessors, and after being notified that a customer is coming, it will try to broadcast on 
channel c immediately with a specific weight oop... The one which succeeds to do so will be the desti- 
nation of the customer and all the others will be informed by receiving a message on c. The other part 
{Wj = 0]0, Cj{wj,Aj) takes care of the requests of the customers in its queue if it is not empty, the rate of 
the customer leaving depends on the service rate. By putting these five servers in parallel, we get the 
whole system CQA/(1,2,3,4,5). Our semantics guarantees that each leaving customer will finally arrive 
at one and only one server Actually, this can be seen as a way to model "point-to-point" communication 
with immediate transitions. For example if a customer is leaving server 3, then both server 4 and 5 will 
be informed by broadcasting a message on channel C3. But after that server 4 and 5 will try to broadcast 
an acknowledge r on channel c, with probability 0.4 and 0.6 respectively. The one which succeeds to 
send r will be the real destination of the customer The other server will know this fact by listening on 
channel c and roll back to its original state at the same time. The model in Table \5\ is quite flexible, for 
example, we can add self loops easily, that is, the destination of a leaving customer can be the same 
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server as its departure server, or instead of having fixed transfer probabilities, we can make them change 
based on the current lengths of queues. 

In Example |3] we have shown a typical closed queueing network where we assume that every server 
has a queue with infinite capacity. But in practice the capacity of queues is often limited. In ifTOll a variant 
of closed queueing network is proposed where the queue of each server / only has finite capacity Bi for 
1 < / < M. A customer which requests service at server / while the queue of server / is full will instantly 
be routed to another server j with probability pij as if it is served by server / at an infinite speed. For this 
kind of model, it is hard (if not impossible) to model it in a compositional way without using immediate 
actions. With immediate actions, the model is easy to obtain without changing the model in Table[5]a lot. 

Example 4 Suppose the queueing network we are about to model is the same as the queueing network in 
Example\3\except that the capacity of every queue is maximum 10. If one server is full, it will just transfer 
the coming customers to its next servers. The model can be obtained by simply replacing the SQj(wj +1) 
in Rev(i,j,Wj) in Table\5\with [w; = 10]c,(ll, ooj), SQy(Wj + 1). This means instead of accepting any 
coming customers, we require an extra checking on the current queue. If a customer arrives at a server 
whose queue is not full, the customer will be accepted, otherwise, the customer will also be accepted 
but will be transferred to the next server via action Ci{ll,ooi)just like it is served with infinite rate, that 
gives excuse of parameter 11. The same process will continue until the customer arrives at some server 
which has free space for it. 

From the semantics in Table |4l we know that all the outputs are non-blocking, i.e., for one message 
to be broadcasted, it is not necessary to have recipients. But sometimes we may have models where 
components are not completely independent and one component can do something only after some other 
components finish, that is, some behaviors are blocking. Blocking here means that an output action can- 
not happen spontaneously but has to wait until certain conditions are fulfilled. For instance in Example 
m every queue has finite capacity. When a customer arrives at a server without free space, it will simply 
be transferred to other servers. What if the customer cannot be transferred but can only wait until the 
server has free spaces? Refer to the following example from lITSl which is also a variant of queueing 
network called open queueing network with blocking. 

The network consists of A'^ parallel servers called merging queues; there is a queue receiving the 
outputs of these merging queues and is called merged queue (or queue 0). The service time at queue / 
is exponentially distributed with rate /I,. The queue network is open since the number of the customers 
circulating in the network is not fixed and some external customers may arrive from the outside. Arrivals 
to queue / are independent Poisson Processes with rate yu,. There is no external amval to the merged 
queue. The length of the /-th merging queue is B,-. The capacity of the queue is So, and Aq is its service 
rate at queue 0. If a customer amves at a merging queue when it is full, the customer will be lost. When 
a customer completes service at server /, it will be transferred to queue only if it is not full; otherwise, 
the customer waits in the /-th server until it can enter queue 0. During this time the /-th server cannot 
serve other customers that might wait in its queue. In this case, the queue is said to be blocked and queue 
is blocking. Since there are M servers in parallel, there might be more than one queue blocked at the 
same time. When more than one queues are blocked, it is assumed that they will enter queue on a 
"First-Blocked-First-Enter" basis. 

From the description of open queue network with blocking, we know there are two kinds of actions 
involving in this model: one is blocking and the other one is non-blocking. For instance, when queue 
is full, any other arrivals have to wait until queue has free space, this is blocking action. On the other 
hand, when the merging queues are full, instead of blocking the external arrivals it will just discard them. 
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Table 6: Model of Open Queueing Network with Blocking 
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so the external arrivals are non-blocking actions. In the following we will show how to model both kinds 
of actions in our calculus. 

Example 5 Fig. \3\gives a concrete example with 3 servers marked as Q[, Q2, and Q3 respectively. The 
length of each queue is 3 and the numbers on the in-edges and out-edges are used to denote arrival rates 
and service rates, that is, Ai = 3, /I2 = 8, A^ = 6, //i = 2, p.2 = 4, and ^t, - 5, and Aq - 10. The Qq at the 
bottom is the queue with service rate 10, and the capacity of its queue is 5. Initially, every queue is 
empty. 

Here is the model of the queueing network in Fig. \3\where Ai denotes the arrival process of queue i, 
Qi{w) denotes queue i with w customers in the queue for 1 < / < 3, and Qo{w,S) denotes queue with w 
customers. The S is a sequence of queues which are waiting for queue when it is full, it is an element 
of SQ which is defined by enumerating all the possible sequences of queues waiting for queue 0. The 
symbol ± is used to denote empty sequence. 

SQ = {±, 1,2,3, 12,21, 13,31, 23,32, 123, 132,213,231, 312,321) 

For simplicity, we define two functions on this set: H{S) : SQ —> {±, 1,2,3} and T{S) : SQ — > SQ which 
return the head of sequence and the left sequence by deleting the first element respectively. For example, 
H(123) = 1 and r(123) - 23; H{±) - ± and T{±) = ±. We use Si to denote a new sequence by attaching 
i to the end of S . Note here that these functions are just used to give a compact model, they can be 
replaced by the standard operators by enumerating all the possible cases. It is similar for conditions like 
{w = 0,1,2 aH(S) = /]. The model of each component in Fig.\3\is shown in Tabled The whole system 
can be denoted as P = \\Ai\\Si{0) \\ 50(0,-1) with 1 < / < 3. 

(' i 

As we said before, broadcasts such that Si{leavej,Ai} are blocking, so when there is no input Si{x,w) 
available, that is, the queue is full, server i has to wait before it can perform other actions. To do so, 
we let the server i wait for the message unblocking on channel bj after one customer leaving from it to 
server 0. If queue has free spaces, it will perform action bi{unblocking,ooj) right after it receives the 
request from server i. Otherwise if the queue is full, it will attach the request to the end of its waiting 
list. When server is ready to handle the request after several steps, it will send the message unblocking 
to server i instantly via immediate action bi{unblocking, 00 j ). Server i will receive it at the same time 
and then unblock itself. 

Fig.^shows a fragment of the execution of P where ^, — ^, and ■>- denote Markovian, 

immediate and probabilistic transition respectively. In additional A =\\i Ai and S{q\,q2,q^) =\\i Siiqi), 
when qi is barred, it means that server i is blocked. For example, S{q\,q2,q-i) = S \{q\)\\ b2(x,l).S2(q2) II 
S^iqs). When in state A \\ 5(3,2,1) || So{4,3), it means that there is a free space in queue while the 
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A||S(3,2,l)||5o(5,3) ^ 



J;: 



A||5(3,3,l)l|5o(5.3) 

■■■ A||S(2,2,l)||5o(5,31) A || S(3,I,I) || So(5,32) 



.■' » \ A||S(3,2,I)||So(4,3) 

y 



♦ i A||S(3,2,l)||So(5,J-) 



Figure 4: Execution Fragment of P 

server 3 is waiting for service. In this case, server should response to it and transfer to the state 
A II S(3,2, 1) II 50(5,-1) instantly where the request is removed from the waiting list to the queue 0. The 
self loop with probability ^ denotes that the arrivals of external customers to server 1 while its queue 
is full. In this case, the arriving customers will be discarded without causing any effects. The other self 
loop with probability ^ is similar except that the arriving customer is discarded because server 3 is 
blocked. 

From Example [51 we can see that blocking actions can be represented easily by using immediate 
transitions. In general when broadcast a{b,A) is blocking, it should be prefixed with an input such as 
c{x,w).a{b,A). When certain conditions ai^e fulfilled, the process should trigger a{b,A) by sending a 
message on channel c instantly, that is, by action Uke c{unblocking, ooj^.). 



5 The Underlying CTMC 

In Examples [3l HI and [51 we see that immediate transition is indeed powerful to model some systems. 
But the main disadvantage is that the underlying CTMC of a process is not so obvious anymore. In this 
section we will show how to define the underlying CTMC in case of immediate transitions. Different 
from 111 where the eliminations of immediate transitions are based on the weak behavioral equivalence, 
we solve this by dealing with a set of equations as follows. 

In this calculus choices between Markovian actions are probabilistic depending on their rates while 
choices between immediate actions are also probabilistic depending on their associated weights. In 
addition, if both types of actions are involved in a choice, the immediate action should be prioritized, 
since they take no time, so the probability of the Markovian action being executed before the immediate 
one is zero iT4|. Since the priority of immediate actions are not shown in Table [H a CTMC cannot be 
obtained directly from a process based on the semantics. In the following we distinguish between MP 
and IP. MP only contains processes which do not have immediate transitions while IP only contains 

processes which have immediate transitions available. Formally, \P = {P e J^ \ 3AV.P - - ^ AP) and 

MP ^ {P e ^ I $A¥.P- - ^AP} ^ ^ \ IP. It is not hard to see that every state in a CTMC should be 
seen as a Markovian process in this calculus, so instead of considering all the processes in Der(P) as in 
Section[3l we only need to consider set Der(P) n MP when talking about the corresponding CTMC of P. 
The question now is how to give the value of M{P, Q) for any P,Q e Der{P) n MP. Due to immediate 
actions, it is not enough to just consider one step Markovian transition as before since P might have to 
go through several immediate processes before reaching Q. 
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Figure 5: Example with Immediate Transitions Figure 6: CTMC corresponding to P 

We use Pr'(R, Q) to denote the probability from Rto Q via all possible immediate transitions where 
Re J^ and Q e MP. The value of Pr'{R, Q) is given by the smallest solution defined by the following set 
of equations: 



Pr\R,Q) = 



{\ R = Q 

Z AR(P)>o AR(P) X Pr'(P, Q) 7? - - ^ AR A Der(/?) n M P ?t (2) 

otherwise 



Then for any P,Qe MP, MiP, Q) can be defined as follows: 

M{P,Q)^Ax Y_^ AP(/?) X Pr'(/?, P— ^AP 

AP(R)>0 

Example 6 Suppose P = a{b, 10) || a(x, l).c{bi,oo2) \\ a{x, \).c{b2,°o-i) \\ c{x, \).c{x,6), by the semantics 

in Table^we can draw a derivation tree as Fig. |6] We omit passive transitions like ^ here. 

In Fig. \5\nondetenninistic choices emerge, such as process c(&2^°°3) II c{b\,6) can choose either 

immediate transition — ^ or Markovian transition ^ . But since the probability of the Markovian 

transition being executed before the immediate one is zero, so the transitions inside the dashed rectangle 
is impossible and should be ignored. 

In this example, Der{P) ^ {P,Pi,P2,P3,P4,P5,Pi,0}, Der{P)r\lP={Pi,P2,Pi,P5,Pi], andDer{P)r\ 
MP = {P,P4,P^,0}. To define the CTMC of P, we only need to consider the processes in Der{P) n MP 
and the corresponding CTMC of P is shown in Fig. |6| which is quite simple compared to the derivation 
tree in Fig. \5\ 

In the second case of Equation (O, we require that Der(/?) n MP t^ 0, that is, there exists at least a 
Markovian process which is reachable from 7?. But sometimes it is also possible for one process reaching 
an immediate state from which no Markovian process can be reached, that is, we have immediate loop. 
We call states in SP - [P \ Der(P) c IP) absorbing states and use a special process Stuck to denote them. 
Accordingly, the set of states of the CTMC should be (DerCP) n MP) U {Stuck} and A^CP, Stuck) where 
P € MP can be defined as follows: 

A4(P,Stuck) = ^- Y_i M{P,Q) P— ^AP. 
QeDer(P)nMP 
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6 Conclusions and Future Works 

In this paper we give a stochastic broadcast k calculus which is useful to model some server-client based 
systems where synchronization is governed by only one participant. Both Markovian transitions and 
immediate transitions are taken into account. A few examples are given to show the expressivity of the 
calculus which is enhanced a lot with immediate transitions. The semantics is given by Labeled Tran- 
sition System without relying on techniques such as multi relations, proved transition systems, unique 
rate names, and so on. Each transition is labeled with rate or weight instead of action and the resulting 
distribution is over pairs of actions and processes instead of only processes. In this way, the underlying 
CTMC can be obtained naturally even with existence of immediate transitions. 

A number of further developments are possible. In the future we would like to provide semantics to 
some of the most representative stochastic process languages as ||9ll4l and compare these different ways. 
Another possible extension is to support parameters, that is, we do not need to know the value of each 
parameter at beginning. We can also reuse the model by assigning parameters with different values and 
so on. Sometimes, some CTMCs have special form, that is, product form which can be solved efficiently 
|[T6l [T2I . We try to answer whether the underlying CTMC of a process is in product form or not by 
syntax-checking. 
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